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y—i 1 Abstract 

o : 

£Nj . We consider the FENE dumbbell polymer model which is the coupling of the incompress- 

ible Navier-Stokes equations with the corresponding Fokker-Planck-Smoluchowski diffusion 
equation. We show global well-posedness in the case of a 2D bounded domain. We assume in 
the general case that the initial velocity is sufficiently small and the initial probability den- 
sity is sufficiently close to the equilibrium solution; moreover an additional condition on the 
(_,' coefficients is imposed. In the corotational case, we only assume that the initial probability 

Oh ■ density is sufficiently close to the equilibrium solution. 
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■ 1 Introduction 

o\ ■ 

t-h . The success of Kirkwood, and of Bird, Curtiss, Armstrong and Hassager (and their collabora- 
Tj- | tors) kinetic theory of macromolecular dynamics triggered a still on-going flurry of activity aimed 
to providing molecular explanati ons for non-Newtonian and viscoelastic flow patterns. This can 



X 



be reckoned from |BAH87] and |Ott06 , for example. The cornerstone is the so called diffusion 
equation, a parabolic-type Fokker-Planck-Smoluchowski partial differential equation, the solution 
of which is the configurational probability distribution function; the later is the key ingredient for 
calculating the stress tensor. 

The simplest polymer chain model of relevance to Bird et al. theory is that of a dumbbell, where 
the beads are interconnected either rigidly or elastically. Although a crude representation of the 
complicated dynamics responsible for the flow viscoelasticity, the now popular Bird and Warner's 
Finitely Extensible Nonlinear Elastic (FENE for the short; see |War72] ) chain model is capable 
in capturing many salient experimentally observable flow patterns of dilute polymer solutions. 
It was therefore quite natural that many researchers took on exploring the fundamentals of this 
relatively simple model (for more on this and related issues see for example [BE94] and |Scli()(i i. 

The aim of this work is to take on studying the momentum-balance (or Navier-Stokes) equa- 
tions together with the constitutive law for the FENE fluid. The latest is obtained by using 
the so-called "diffusion equation" , practically a Fokker-Planck PDE, the solution of which is the 
configurational probability density. Put it differently, we focus on a system of equations that 
consists of a "macroscopical" motion PDE and a "microscopical" Fokker-Plank-Smoluchowski 
(probability diffusion) PDE. More precisely, given a smooth bounded connected open set ftcK d 
and some ball D(0, R) we will study the initial boundary value problem which consists in finding 
u = u(t, x) : R + x tt ->• R d , g = g(t, x, q) : R+ x Q x £>(0, R) -> R and p = p(t, x) : R + x Q ->■ R 
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solutions of the two following coupled equations: 



7 a.. , n. 7(1-7)^ (f 9®J7 

ill 
R 

and 



1.1) d t u + u-Vu-^-Au+Vp= ' v './ Vs- / -^—^9(t,x,q)dq\ on R + x Q 

Re ReWe^ \Jd(o,r) 1 



1.2) % + « • V x g + V ¥ ■ (a(n)^) = — £^ A 90 + ^e V « ( ^ZW 9 I on M + x fi x ^( ' «)■ 
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Moreover, the vector field u must be divergence free and g must be a probability density in the q 
variable: 



(1.3) diy r w = 0, / gdq = 1, g > 0. 

Jd(o,r) 

The boundary conditions are 

(1.4) M |^= 

plus some boundary conditions for g on fi x dD(0,R) which will be embedded in the function 
spaces we will work with. 

The constant 7 belongs to (0, 1), Re and We are (respectively) the Reynolds and Weissenberg 
numbers and N, R are some polymer related physical constants used to obtain dimensionless 
quantities. We assume all these constants to be strictly positive and moreover that NR 2 > 2. 
The quantity a(u) is a short-hand notation for either Vu or Vu — (Vu)*. In fact, the physical 
significance is achieved when cr(u) = Vu; we will call this the general case. The choice cr(u) = 
Vu — (Vu)* is very close to being physical significant while having better mathematical properties; 
we will call this the corotational case. Let 

1 and - * 



) Id(o,r) ^ 

It is not hard to observe that the couple (0, Z) is a steady solution of fll.ip - fll.4p . 

The initial boundary value problem ( ll.ip -( fO|) was studied by several authors but mostly in 
the case where Q = M? or M 3 . The results are different, depending on the model (general or 
corotational). We start by describing the results where Q = M. 2 or R 3 . We restrict ourselves to 
the model described above, but we would like to mention that there are other results on closely 
related problems (for example a model when the variable q lies in the full plane or full space, the 
Hookean model, etc.). We refer to |MaslO] for a discussion of all these models. 

Global existence and uniqueness of strong solutions of problem fll.ip — (jl.4p is known in the 
following situations: 

• fi = M 2 and corotational model if u e H S (R 2 ) and g G H S (R 2 ; H%(D(0, R))), s > 2 (see 
[LZZ08J). The regularity of go in the q variable was improved in |Mas08j to some LP weighted 
space for large p. 

• Q = IR 2 and general model or Q = M 3 and general or corotational model if Uq is small in 
H S (M. 2 ) and if \\Z~2 \\g — Z\\ H sr R 2\\\ 2 is small, where s > 1 + | where d G {2,3} is 



\L 2 (D {0,R)) ° ^""5 ""^^ " x 1 2 
the space dimension (see |LZ08t IMas08j ). 



Global existence (no uniqueness yet) of some weak solutions for rough and arbitrarily large 
initial data was proved in both dimension 2 and 3, first in the corotational case by |LM07] and 
quite recently in the general case by [MaslOj . see also |BS10j for a slightly different version of the 
system of equations. 

Long time asymptotics of the general model were studied in |JLLO06] where a priori esti- 
mates are obtained to prove formally the stability of the equilibrium solution. In [BS SQ5] the 
authors studied a related model where a smoothing operator is acting on the velocity field and 
the corresponding stress tensor. 

As far as strong solutions on domains with boundaries are concerned, we are aware of two 
works. One is |ZZ06] where local existence and uniqueness is proved if u G H 4 (Q) and if g is H A 
in x and has some weighted if 3 regularity in the variable q. Another one is [KP10] where local 
existence is shown if uq G W 1,p (Q) and g^ is W 1,p in x and has some weighted LP regularity in the 
q variable and p > d. 

The goal of this paper is to address the issue of existence and uniqueness of strong solutions 
for the above mentioned initial boundary value problem on bounded domains Q C M 2 with ho- 
mogeneous Dirichlet boundary conditions. This is not a straightforward adaptation of the known 
results in the full plane. Indeed, the proof of global existence results of solutions proved by [LZZ08J 
uses heavily the Littlewood-Paley decomposition and paradifferential calculus; this is of course not 
available on bounded domains. Even the global existence results for small data involve technical 
difficulties that make necessary to assume an additional condition of the material coefficients, 
more precisely we will need to assume ( 11. 5p . We refer to Section [6] for a detailed explanation why 
this is necessary. 

In the general case, we show the following global existence and uniqueness result for initial 
data which is sufficiently close to the equilibrium solution (0, Z). 

Theorem 1 (general case). Lets G (l,f). A ssume that uq is divergence free, vanishes on dVt and 
belongs to H S (Q). Assume moreover that Z~* ||^o||fr s (n) £ L 2 (D(0, R)) , g > and J D ^ Q R ^ g dq = 
1. There exists two positive constants K\ = Ki(Q, s) and K 2 = K 2 (fl, s, 7, Re, We, N, R) such that 
if the fluid related coefficients verify the relation 

(1.5) < Kl 

V ; NWe ~ 

and if the initial data is sufficiently close to the equilibrium solution (0, Z) 



|«o||fp(n) < K 2 and 



\9o ~ Z\\ H '{q) 



'Z 

then there exists a unique solution to system (ll.ll) - (ll.4p such that 

u G L°°(R + ; H S (Q)) n L 2 (R + ; H S+1 (Q)) 

and 



< K 2 

L 2 (D(0,R) 



> 



Mk s (n)|| L2(D(0jK) G £°°(K+), \\Z*\\Vq{g/Z)\\ H °(n)\\ L2imR) G L 



In the corotational case we improve the previous result in the following manner. Not only the 
restriction on the material coefficients ( II. 5p is no longer required, but the initial velocity Uq is 
arbitrarily large as well. More precisely, we have the following theorem. 

Theorem 2 (corotational case). Let s G (1, |). Assume that Uq is divergence free, vanishes 
on dVt and belongs to H s (Vl). Assume moreover that llgollf^^/V^ G L 2 (D(0, R), go > and 



Jnfofi)^ = ^ ere e2; 2s£s a positive constant K 3 = K 3 (Q, s, 7, Re, We, N, R) such that if the 
following smallness assumption holds true 



\9o 



\H s (n) 



< 



L 2 (D(0,R) 



exp 



K 3 (l + \\u \\ Hsm )e K3llu <^ 



then there exists a unique solution to system fll.ll) - fll.4p such that 

u e L°°(R + ; H S {Q)) n L 2 (R + ; H S+1 (Q)) 



and 



\Z 2 \\g\\H"(n)\\ L 2 



(D(0,R) 



\zv\\vmz)\\ h . 



( n )\\L 2 (D(0,R) 



Compared to the result of |LZZ08j valid in the case of the full plane, we have an additional 
condition on g : it needs to be close to Z. As explained above, this is due to the fact that we 
work with bounded domains and the methods of |LZZ08j do not work here. Nevertheless, we have 
an improvement in the regularity assumptions. More precisely, we require a regularity in the x 
variable which is H s , 1 < s < |, while in [LZZ08| IMas08] it is necessary to assume that s > 2. 
The regularity in the q variable is also improved, roughly from H 1 to L 2 . 

The paper is organized as follows. In Section [2] we reformulate the problem f ll.ip - fll.4p and 
introduce the notations. We construct next in Section |3] a sequence of approximate solutions. The 
global existence of the approximate solutions is proved in Section HJ We show uniform estimates 
for the approximate solutions and complete the proofs of Theorems [1] and |2] in Section |5j The last 
section contains two final remarks on the hypothesis we have to assume. 



2 Notations and functional framework 

We start by making a change of functions allowing to rewrite the equations in a better form. 
Notice first that 



If we set 



and 



V^g + N—^ r g = ZW 1I ( | ). 



1 \1? 

1 R 2 
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q=± M{q) = {l-\q 



2\6 



f(t,x,q) = g(t,x,Rq) 



j_ 7(1-7) (25_ 

ai Re' 0(2 Re We 2 \N 



1 



NR 2 
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then the couple (u, /) must verify the following system of equations: 

q® q 



(2.1) 
(2.2) 



d t u + u ■ Vu — aiAu + Vp = a2V x 



.-/ dg on R + x Q 



d t f + u- V x f - a 3 V q 



MV„ 



f_ 

M 



+ V 9 ■ (a(u)qf) = 

onR + xOx D(0, 1) 
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and 

divu = 0, / fdq = 

Jd{o,i) n 

The boundary conditions are homogeneous Dirichlet conditions for the velocity u plus some bound- 
ary conditions for / on Q x dD(0, 1) which are implicit from the condition that / has the below 
defined H\j regularity in the q variable (we refer to [Mas08j for a discussion on the boundary con- 
ditions verified by functions in H]^). We also prescribe the initial data u\ t=Q = u and / = f . 

We recall now the usual spaces for the solutions of the Navier-Stokes equations: 

H = {v G (L 2 (n)) 2 ; V-v = 0, v ■ is = on dVl} 

where v is the outward normal to dQ and 

V = {ve (H*(n)f; V-v = 0}. 

We will abbreviate in the following D = D(0, 1). We introduce next the following Banach spaces: 

i r 2 

L 2 M = L 2 (D,-dq) = {tp; J^dq<oo}, 

L 2 X L 2 M = L 2 (n x D^dxdq) = {(p; ff ^dxdg<oo}, 

J J Q x D 

L 2 x L 2 M = L 2 (R 2 xD,±-dxdq) = {ip; [[ ^dxdg<oo}, 

M J JR^xD M 

H x L li = W : x D -> c ; IMIff^ ■= \\\\v\\H°(n)\\ L 2 M(D) < oo}, 
H l L li = {V ■ x D C; IMb-z^ := || llvlk'(ii a )|| L a tf( £,) < °°}. 
H° X H X M = { V : n x D -> C; |M|^ == || ll MV ?(^) ll*"(n) Lyi>) < °°}' 
#*#M = : x D C; IMI^ := ||||A^V ff (^)llH-(M>)||^ (D) < oo}, 

H x H M = H x L M n H x H M- 

The quantities (function spaces, vector fields, etc.) with a bar on top have the x variable in R 2 . 

Below, all functions / and their different versions (f n , f n , g n , g n , etc.) are assumed to belong 
to spaces which are in the q variable (which implies boundary conditions in the q variable). 

We have the following Poincare type inequality. If <p = (p(x, q) is such that j D ipdq = 0, then 

(2-3) hhiLi < cy\\ L 2Hi M , 

see jChulOj . 

Let now denote by A£ the Fourier multiplier (D) cr , i.e. the operator of multiplication in the 
Fourier space by (1 + |£| 2 )^. We will always apply this operator in the x variable (the functions 
need to be defined on M. 2 of course). 

The divergence of a matrix is taken along rows: for A = (a^) we define div A = • dja^i. 
Scalar product of matrices is defined by A : B = a ij^ij- The tensor product of two vectors x 
and y is the matrix x ® y = (xiyj)ij. 



3 Construction of a sequence of approximate solutions 



Clearly C °° (D) is dense in L 2 M and a positive function in L 2 M can be approached by a sequence of 
positive smooth functions (cut-off and convolution preserve the sign if the cut-off and convolution 
functions are non- negative). Moreover, we can assume that the integral of each of the approximate 
functions is equal to the integral of the limit function (otherwise we can normalize each of the 
approximate function by multiplying with an appropriate constant). In the x variable, the usual 
standard smoothing procedure is also achieved by cut-off and convolution and these operations 
preserve the sign. Therefore, there exists f$ G C^°(Q x D), non-negative, such that 

(3.1) fZ^fo 
and 



in H S X L 2 M 



(3.2) 



f f£dq = R- 

J D 



Let us denote by P the Leray projector, i.e. the orthogonal projection operator from (L 2 (Q)) d 
onto H, and by A the Stokes operator defined by A = —PA. It is well known that A is a self- 
adjoint operator on H with compact inverse and that D(A a ' 2 ) = (H a (tt)) 2 n V, Ver G [1, |) and 
D(A a / 2 ) = (H a (n)) d H H, Vcr G [0, |) with equivalent norms (see |FM70] ). 



Denote Ai, A2, . . . , A n , . . . the sequence of eigenvalues of A and v\, V2, ■ ■ ■ , v r , 



the corre- 



sponding eigenvectors that form an orthonormal basis in H. Let H n := Jz?{vx,V2, ■ ■ ■ ,v n } be the 
vector space spanned by the first n eigenvectors of A, and P n the orthogonal projection of L 2 (Q) 
onto H n . We endow H n with the L 2 norm making it a Hilbert space. We observe that for any 
a G [0, ~) there exists some constant C(a,fl) independent of n such that for any g G H a (Q) we 
have 

(3.3) \Kg\\ H °(n)< c{a,n)\\ g \\ H , m . 

Indeed, it is well-known that P is bounded on H a (tt). Then P# G H a (Q) n H = D(A a l 2 ). 
Moreover, P n is an orthogonal projection in D(A°~/ 2 ) so we can write the following sequence of 
estimates: 



Kg\\H*(n) < C\\A°/ 2 F n Fg\\ LHn) = C||P n A CT / 2 P^|| i2 



»(n) 



< C\\A°/ 2 Fg\\ LHQ) < C'\\¥g\\ H , {u) < C"\\g\\ H 



which proves ( 13. 3p . 

Letting P operate on relation (12. ip leads to 



du 

— + a x Au 
at 



[u ■ Vu) = a 2 P 



D 



1 



\Q\ 



fdq 



that 
(3.4) 

and 
(3.5) 



We consider the following approximation problem: find (u n ,f n ) with u n G C([0, T], H n ), such 

q® q 



d t u n + ai Au n + ¥ n (u n ■ Vu n ) = a 2 F n 



D 



i-M 



;f n dq 



d t f n + u n -V x f n -a 3 V, 



MVr, 



r 

M 



+ V g -(a(u n )qP) = 



with respect to the initial conditions 
(3.6) u r 



\t=Q 



ul = F n u Q , 



r 



\t=0 



and such that /" has H\j regularity in the q variable. We will later use that 
(3.7) < u in H s (n). 



4 Global existence of the approximate solutions 

Let us first remark that for any / G H n and any m G N, one has 
(4.1) ||/|k™ { Q)<C(^m,n)||/|U 2(n) . 
Indeed, if / = Y17=i aiVi then 

/ n \ 1/2 

fl» (n) < C||A^/ 2 /IU 2( o) = C $>, 2 A™ < C u*x\? /2 \\f\\» (n) . 



,i=i 



Throughout this section C denotes a constant that depends on n, Q, material coefficients and 
other constants. It may change from one line to another. 

The existence of the approximate solutions is granted by the following theorem: 

Theorem 3. There exists a global solution (u n ,f n ) to the problem (13.41) - (13.61) . such that u n G 
C°(R + ;H n ). 

Proof. In this proof the various constants C, Ci, C2, ■ ■ ■ may depend on n but are independent 
of time. We fix an arbitrary finite time T > and show that we can solve (I3.4l) - (l3.6p up to time 
T such that u n G C°([0, T\; H n ). 

Suppose that u n is an element of C°([0, T]; H n ). From (14.11) we deduce that 

sup < +00. 

*e[o,r] 

Therefore one can construct the flow Xn{t, y) of u n as the unique solution to the equation 

d t Xn(t, y) = u n (t, Xn(t, y)), x«(0, y) = y. 

Since the Jacobian determinant of Xn(t, •) is equal to 1 we deduce that for any t G [0, T], Xn(t, ■) 
is a C°° - diffeomorphism from Q to Q. 

Let f n (t, y, q) = f n (t, Xn (t, y),q). Clearly f n solves (133 j) - ([3TB]) if and only if f n solves 



(4.2) 



'd t f n + V q ■ [a(u n ) o Xn qf n ] ~ a 3 V q ■ [MV q (| 

f n (o,y,q) = K(y,g). 







In the equation above, the variable y plays the role of a parameter only. The existence, 
uniqueness and smoothness of a solution f n to (14. 2p which is H\ { in q was proved in [Mas08j . This 
allows to construct f n , and therefore f n , if u n is given. We denote by S the operator that gives 
f n in terms of u n , f n = S(u n ). 

Since u n is smooth enough w.r.t. x we also deduce that f n is smooth enough w.r.t. y and the 
same holds true for f n . Then f n satisfies (I3.5P and (13.61) . We observe moreover that 

(4-3) // n d^4- 



D 



R 2 



Indeed, from (13.21) we know that the above relation is satisfied at time t — 0. If we integrate with 
respect to q relation (13. 5ft we have that the quantity J D f n dq is transported by the vector field u n 
so it must be constant. 

We conclude from the preceding observations that it suffices to show that there exists a global 
solution u n G C°([0, T]; H n ) of the following equation 



(4.4) d t u n + F n (u n - Vu n ) + ai Au n = a 2 F n (v x ■ j S(u n )F{ 



(q)dq), u n (0,x)=F n u 



where F(q) is the following matrix: 

F[q) = T^W' 

Indeed, if such an w n is obtained, then the couple (u n , S(u n )) solves (13.4p - (l3.6p . 

To solve f)4.4p we will use a fixed point method. More precisely, we write ( 14. 4 p under the 
following equivalent integral form: 

(4.5) u n (t) = e- aitA F n u + [ e ai{s ~ t)A F n [-u n ■ Vu 11 + a 2 V x ■ [ S(u n )F(q) dq] (s) ds 

J0 JD 

We search for u n as a fixed point of the operator 

B : C° ( [0, T ] ; H n ) C°([0, T ] ; H n ), B(v) is given by the rhs of (Q 

where the time T is to be chosen to ensure B is a contraction mapping. Recall that H n is endowed 
with the L 2 topology. One has: 



(4.6) \\B(v)-B(v') 



\L 2 (n) 



< 



e ai(s - t)A F n (v ■ Vf -v' ■ Vv') 



L 2 (Q) 



+ « 2 



e a ^ A F n [V x - f (S(v)-S(v'))F(q)dq] 
Jd 



h + h- 



L 2 (n) 



Next: 



h 



< 



(4.7) 



/)-VuW-V(?;-?/)|| L2(n) 



L 2 (Q) 



f — f 



< C / ||u - T/||fla(n)(IM|fla(n) + lk'lli? 2 (n)) 

Jo 

< Ct SUp ||U - t/|| L 2 (n) (sup |M|L2 (n) + sup ||u'||L2(n)) 

[0,t] [0,t] [0,t] 

where we used relation ( 14.ip . To bound we observe first that we have the following inequality 

n n n 

l|P»V • ^|| 2 L2(n) = E K V • h ^i)\ 2 = E K*> V ^)! 2 ^ E ll^ll^njIlV^H^n) < C(n) 

t=0 

Therefore 

J 2 < «2 



2 

L 2 (Q)- 



i=Q 



8=0 



P„[V«b- / (S(v)-S(v'))F(q)dq\ 
Jd 



L 2 {Q) 



(4- 



< Ctsup 

[0,t] 

< Ctsup 

[0,i] 



(S(«)-S(«'))F( ff )d« 
S(v) - S(v' 



D 



L 2 (nxD) 



L 2 (p.) 

MF\\ L 2 {D) 



<Ctsu V \\S{v)-S{v')\\ LlLli 

M 



s 



where we used that S > 1 to have || "vMF\\ L 2^ < oo. Next, we remark that $ = S(v) — S(v') 
solves the equation 

9 t $ + (v- v') ■ V x S(v) + v 1 ■ + V g • [a(v - v')qS(v)] 

+ V g • (a(v')q®) - a 3 V q 



M 



0, $(0,x,q) = 0. 



We multiply the above relationship by — and integrate in x and q to obtain, after some straight- 
forward calculations that 

d „ , „ „ „ T „o „/"/",, , , $ 



+ 2a 3 \\n 2 LlHli = 2 1 1 {v'-v)- V x S(v)- + 2 / / a(v - v> 



+ 2 /Mt;'): 



M 



17 



17 



<2||^ — v'lU-c^llv^s^)!!^^!!^!!^!^ + 2||cr(t; — ^011^(^)118(^)11^^11^11^^ 

+ 2|k(t;')|U-(n)||<f||L S Ll / ll < f|lLli/l f - 

As in (14.31) we have that j D $ dq = so the Poincare inequality (12.31) holds true for We deduce, 
using inequality (14. ip . that: 



(4.9) -\\n 2 L2L 2 +a 3 \\m 2A1 <C\\v-r 



/ 1 1 2 



^(iD [II V*S(<,)||^ 2m + ||S(t;)||^J 

+ ^IMIiWI $ llW 



M 



0, S( W )(0,x,g) = /™(x,g) 



Recall that S(v) solves the equation 
(4.10) d t S(v)+vV x S(v)+V q -[a(v)qS(v)}-a 3 V q - MW q 

S (v} 

Multiplying the equation of S(v) by — and integrating w.r.t. x, q, we obtain after similar 
estimates the following relation: 

d 



(4.11) 



(l J|SWII^ + « 3 ||S( U )||^<C'||,;||i 2(n) ||S^)||i, i2M 



We will use in what follows several times the following simplified version of the Gronwall 
inequality: 

y' < aiy + a 2 , y(0) = y 
with ai,a 2 positive constants, implies 



y(t) < (y + a 2 t)e a ^. 



From (Kill it follows: 



(4.12) | 
where we denoted, for any z G 



SM(t)\\ 2 LlLlt <Rl\ sup |Ml!a ( n) 

tG[0,T] 



= ll/olliiz^exp (CzT). 



In -^ m d .S(v) 

Similarly, differentiating (I4.10p with respect to Xj, multiplying by — and summing over i, 

gives after integration and some straightforward estimates that 

^l|V,S(,)||^ Li/ + 2a3||V,S(,)||^ 

<C1|Vv|U»(n)||V a S(«)||2| +C||a(i;)|| L o 0( n)||V«S(i;)||^||V x S(T;)|| La ^i 



C r ||V<7(v)|| L -(n)||S(t;)||^ £ a f ||V a .S(«)|| i a A i f 



<a 3 \\V x S(v)\\ L2 
Using (14.121) and the Gronwall lemma implies that 



(4.13) \\V x S(v)(t)\\l 2Lh <R 2 2 \ sup \\v\\ 2 L2m 

te[o,T] 



where 

r2 r2 



= [l|VJolI |LL + C*i2?(*)Tj exp (C(z + 1)T). 
From (Oil . (Q2]) and (l4~T3D we get 

l|S(«) - 8(10115^ < CT||t; - vTw T [#L (Nik) + R\ [Ww T )\ **P {CT\W\\w T ) 
where we denoted 

W T = C([0,T];H n ). 
Combining the above with (14. 6p . (14. 7p and (14.81) yields 

(4.14) \\B(v) - B(v')\\ Wt <C x T\\v-v'\\ w T i\\v\\w T + \\ v \\w T ) 

+ ai*' 2 \\v - v'\\ Wt [Rx(\\v\\ 2 Wt ) + MPWwr)] ex P ( C AW\\w T T) 
By a procedure similar in nature to the one detailed right above, one gets: 
(4-15) \\B{v)\\ Wt < IKH^+^Tl^ll^ + CaT^dl^ll^) 

Let now Kq be such that 

max (K||i2(n), H/olL^ , W^xfohiq,) < K . 

Now taking a fixed sq such that so > K (for example so = -Ko + 1) we can choose a T sufficiently 
small such that 

Ko + CzTosI + CzTqR^sD^so 

and 

2CiT s + Ci^o 72 [Ri(sl) + R 2 (sD] exp(C 1 ^T ) < 1 

The above assumptions together with (14.141) and (14.151) show that the operator B is a contraction 
from the closed ball -8(0, sq) of Wt onto itself. The fixed point theorem can therefore be applied 
to grant the existence of a local in time solution on [0, To]. Moreover, the local time existence To 
depends only on the bound K for the initial data Uq and S(t>)(0) = fft. Starting from time T , 
the same argument can be applied to extend the solution, and so on. We justify now that T can 
be reached in this way in a finite number of steps. When re-applying the fixed point argument 
from time T , the new time of existence depends only on max (||ii n (T )||£2(Q), ||/™(T ) U/ji^a ). But 
the estimates shown in the next section imply that any solution (u n (t), f n (t)) on a time interval 
included in [0,T] can be bounded in L 2 x H\L 2 M independently of t G [0, T]. Indeed, we will show 
in particular some H S (Q) bounds on u 11 . This implies L 2 (Q) bounds on u n and by estimates (I4.12p 
and (14.131) some H],L 2 M bounds on f n . This means that the time-existence To can be chosen the 
same at each step, so the time T will be reached in a finite number of steps. This completes the 
proof of Theorem [3J □ 



in 



5 Uniform estimates for the sequence of approximate so- 
lutions and end of the proof 



Let us introduce the new function 



r = f n - aM{q) 



where 
(5.1) 



5 + 1 



R 2 f D M(q) dq nR 2 



It is not hard to check that the couple (u n , f n ) verifies fl3.4p - fl3.6l) if and only if the couple 
(u n , ip n ) verifies the system of equations (I5.2l) - ()5.4p below: 



d t u n + ai Au n + F n (u n ■ Vu n ) = a 2 P r , 



D 



q®q 
1 — |gp 



ij n dq 



fip T 



mv 4m 



+ V, • (a(u n )qi} n ) = -aV g • [a(u n )qM] 



(5.2) 
and 

(5.3) d t i> n + u n ■ V x ^ n - a 3 V g • 
with respect to the initial conditions 

(5-4) u n \ t= = u n , r\ t=0 = W = fo ~ aM(q). 

We observe that relation (13. II) can be rewritten as 

(5.5) Vo^^ Q = /o-aM in H S X L 2 M 
as n — )• oo, while relation (14 . 3 [) is equivalent to 

(5.6) / i) n dq = Q. 

J D 

Remark 4. We observe that 

V q ■ (a(u n )qM) = V q ■ (a(u n )q)M + (o-(u n )q) ■ V q M = (a(u n )q) ■ V q M. 

We used above that tr[a(u n )] = so V g ■ (a(u n )q) = 0. Since V ? M is proportional to q we infer 
that if a(u n ) is skew- symmetric then 

V q ■ (a(u n )qM) = 0. 

We conclude that in the corotational case the right-hand side of the equation ofip n given in relation 
(15.31) vanishes. 

We will show now some uniform (in n) estimates on the approximate solutions u n and ip n 
constructed above. In the calculations below, C is a generic notation for a constant that does not 
depend on n and its numerical value changes from one calculation to another. 



1 1 



L 2 estimates in the corotational case. Our first bound is a L? energy estimate on 



Lemma 5. We have that 

K(*)lli»(n)+«i/ \\^\Un)<\\MUn) + ^4Ml^ Vt > 0, 

where 

- TO 2 

85 4 aia3 

Proof. This estimate is well-known in the case of domains without boundaries. We only have to 
check that it goes through in the case of domains with boundaries, and also that it is compatible 
with our approximation procedure. Moreover, we need to compute the precise best constants so 
we will give a detailed proof. 

Recall that a(u n ) = Vu n — (Vw n )'. Given that a(u n ) is trace- free we have that 

V q ■ {a(u n )q^ n ) = a{u n Y : q <g> V^". 

We multiply next (15. 3p by tp n /M and integrate in x and q. After recalling that the right-hand 
side of (15. 3p vanishes in the corotational case, we obtain 



ld ,L,n„2 , „. , L ,„„2 _ ff -^*. g ^ V ^n# 

a(u n f :q®V q {\4> n \ 2 )M 



1 

(5.7) "2 

1 



2 
0. 



cr(ii n )* : V,M|^ n | 2 

S!xD 



since cr(tt n ) is trace free and skew-symmetric and g £g> V q M is symmetric. 
We multiply now (15. 2 j) by w n and integrate in space to obtain 



1 d 

2dt"" lli2 (°) 



"™ll« n ll? 2 ^ + «i||V M n ||l 2(n) = a 2 / u^V*- / -i^^dgda; 

Jn Jd - 1 - — |9| 



-a 2 // V«»:^" 



q®q , n a 2 ff _ n _ . . ^ n a 2 



Vu n : g$ 






■ ™ l lb 




8<5 4 ai 11 





(5.8) 



Next, 

< ^l|V«"||y (I »||?Vji?|U.( D) ||V."|| l|Si( < yl|V«"|ll I(n) + 
We used above that 

(5.9) \\qVM\\ LHm = -. ^ < 

We infer now from ( 15. 8 p that 

sii«ii. m +«.iiv«ii, m) <^ r ii*"iiv„- 
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Integrating the above relation and (15. 7p in time implies that 



« B (t)i& (0 ,+ai/ nwii 2 L2(n) < H\\i 2{n) + ^-ur 2 







This completes the proof of the lemma. □ 
Remark 6. A similar L 2 energy estimate holds true in the general case too, see JJLLO 06J. 

H s estimates for u n . The next step is to prove H s estimates for u n . 
Lemma 7. There exist a constant C\ > depending only on Q and s such that: 

(5.10) 
d 

df 



4l|A f M"||| 2(Q) +a 1 \\A 1 t s u n \\l 2(n) < C 1 min[/i(0||Aiw n ||| 2(n) ,||A^ n || Z 2 (n) ||A^« n ||| 2(n) ] 



+ C l ^\\A^u n \\ L , m \\r\\ Hg H li 
where the function h is integrable on R + and satisfies 

(5.11) / ^)^<«i 7 (ll^oiii 2 (n) + «4|i^oiii iL |J fl . 

J o 

Proof. Let us take the scalar product in L 2 (Q) of ( 15.2ft with A s w n . We get 
^l|Ait» B ||i a(n) + a 1 ||A^t»l^ (n) = - / A^P n (n"-V^)-A^n"dx 

Jn Jd 1 — |9| 

Using that s G (1, |) and relation (13. 3 j) we infer that 

(5.12) ^||Ai^||| 2(n) + 2a 1 ||A i *V||i 2(n) < C|K • W1^-x || A^ U "|| L2(n) 



+ a 2 ||A^P4V,-^ I ^^dg] 



II A nil 

| L 2 (n) l|A 2 U ||i2 (n) . 



Standard product rules and interpolation in Sobolev spaces imply that 



1 S_ S_ I £ 

IL, n Y7„, n ll r^\\n. n \\ 1 1 _ .n 1 1 ^ ✓"'IL.fMI 2|L,n||2 II,., nil 2 IL,«II 2 

||w ■ v« ||hs-i<c-||u ll//§H n IIh 1+ § — MI u II i2 |p IIhiIp ll#«|p llifi+ s 

We infer that 

|1 — ^11 rj ll ^ II a — nu^ II a iil «ii2 — 



2 



C|K- VM n || H -i||A— M n || L2(n) < CII^II^^II^II^IIA^II^HA— «1z/ 

< axllA^^llla + Cal'-IKIllr^lV^lliallAVllia. 
On the other hand, using again the product rules in Sobolev spaces we can also estimate 

\\u n ■ Vw n ||^-i < C\\u n \\ H s\\Vu n \\ H s-i < C\\u n \\ 2 Hs < C\\A-2u n \\ 2 L2 < C\\Aiu n \\ L 2\\A^u n \\ L 2. 
Hence 

(5.13) C||w n • Vu n || H .-i||A 1 * £ u n || La( n) < a^A 1 ^ u n \\ 2 L2 

+ Cmin^^^Hw^lll^ 2 !! V^Hl^lA^^Hla, ||A5 W -||^ 2 1| || | 2 ] 
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To bound the last term in (1 5 . 1 2 [) we observe that 

V ' ' L T^? r Aq So ^ K * (?M)1 d? + ^' ■ [H So Ti d «l 
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where Id denotes the identity matrix. The last term above is a gradient, so it belongs to the kernel 
of P„. Therefore, the last term in (I5.12p may be estimated as follows 



||A^P„[V.. J D ^^rdq]\\ L2 ^ = l 5 \\A^F n [V x - jv q Q® q M& q ] 

(5.i4) <y!iril^!kVM|U 2p) 

where we used (15.91) . 

Relation flBTTO]) follows from relations ( l5TT2j) . fl5TT5|) and (01)1 if we set 



L 2 (Q) 



h{t) = Ca\ '||M n ||| 2 2 ||VM n ||2 2 

for some suitable constant C . Relation (15. lip is a consequence of Lemma and this completes 
the proof of the lemma. □ 



H s estimates for lp n . We need now estimates on ip n . For technical reasons, in order to 
obtain these estimates we need to work in M 2 for the x variable. Because there are no boundary 
conditions for if) in the x variable, it is possible to extend the equation of tp n to M 2 x D . 

Let E be a total extension operator from Q to M 2 , i.e. a linear operator bounded from H a (Q) 
to H a (M. 2 ) for every o > (in fact we only need it to be bounded for a G [0, 3], that is we only 
need a 3-extension operator). The existence of such operators is well-known, see e.g. |Ada75| 
Chapter 4]. 

We define now 

^ = Em 

so that tpQ G H S X L 2 M and 

(5-15) lKlk^ M < C o\M\h^ m , 

where Cq depends only on Q. The way the extension operator is constructed in [Ada75] ensures 
that the integral in the q variable is preserved: 

(5.16) / V^dg^O 

Jd 

since (I5.6P holds true. 

Next, we want to extend u n to a smooth divergence free vector field defined on IR 2 . In order to 
preserve the divergence free condition, we need to introduce the stream function. For a divergence 
free vector field v defined on Q and vanishing on the boundary of Q, it is well-known that there 
exists a stream function, i.e. a scalar function J such that v = V^J. Moreover, since v vanishes 
on dfl we have that J is constant on each connected component of dQ. Let To be such a connected 
component. Since Q is connected, we clearly have existence and uniqueness of J if we impose that 
J vanishes on Tq. In the sequel, we define 3(v) as the unique stream function of v vanishing on 
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IV By the Poincare inequality, we have that J is bounded from H a (Q) D Hq(Q) to H a+1 (Q) for 
all o > 1. 

We define now W 1 = \7 ± E(J(u n )). Clearly M n | Q = u n and 

(5.17) \\u n \\ H . m < C\\E(J(u n ))\\ H , +HM2) < C\\3(u n )\\ H ^ {n) < C\\u n \\ H '<ti) Va > 1, 

where C = C(Q,o~). 

We finally define ip as the unique solution in R 2 x D of the PDE 



(5.18) d t T + u n - V x T-a 3 V q 



MV ^M 



+ V g • (a(u n )qT) = ~oS7 q ■ [a(u n )qM] . 



The existence and uniqueness of such a tp follows from the argument given at the beginning of 
the proof of Theorem [3] (the variable x plays the role of a parameter only) . By uniqueness of 
solutions of (I5.3P and 05 .4|) . we have that ip \ nxD = ip n - Moreover, given (15. 16[) we can prove as 

for ip n that relation (15. 6 P holds true for ijj n \ 

(5.19) / ^dg = 0. 

Jd 

The following lemma gives our estimates on ^T. 
Lemma 8. There exist a constant Ci > depending only on Q and s such that 

(5-20) ^lirilW + 2a s ||rilW < ^IIA^u-ll^liril^lirilB-A 



m t/ie general case and 



— 1| A 2 u H^njll^ 



(5-21) A||r||^+2a 3 ||? l ||^ < ^||A^|U 2(n) |rik^||? l |k^- 

m t/ie corotational case. 

Proof. We apply the operator A* to (I5.18p . multiply by A x ip n /M and integrate in x and g to 
obtain 



J JR 2 xD ^ 

■ / / A^V, • [a(n") ff M] 
J iK 2 xfl 



M 



= A + J 2 + J 3 



We bound first /j. Let [A*, if 1 ] be the standard commutator defined by [A s x ,u n ]f = A^tT 1 /) — 
vPK. s x f . Using that it™ is divergence free, we can write 



« 2 xD J JR 2 xD 1V1 J JR 2 xD 

< ||A^|| r2i2 ||div.[Ai,B»]^|| 2?£SI < \\KTy L2 \\\\[K,u n }T\\m 



<C\k x ^) || r 2 L2 ||||'U n || H3 +l( R 2)||'?/' ||h»(R2)|| l 2 = CII^Hho+^K 2 )!!^ Ilff 



2 



<C||A^n«|| L2( n)||ril|^ 2 
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where we used the embedding H s C L°°, the classical commutator estimates, see |Tay91 , Section 
3.6], and relation (I5.17p . 

Next, we write 



T 2 

M 



< CW^WhIh^W 11^(^)11^(^)11^11^(^)11^ < C , ||^ l ||^^||^r||^a L ^||A 1 t s w n || L 2 (n) . 



We now make an integration by parts in I 3 and bound as follows 
h = a !! A:[a(t^) ? M]-V 9 (^) <a||AM^)IU W ||?VM|| i2 

JJm. 2 xD M 



<^I|A» M ' 



U 2 (n)ll^ \\HtH\ 



* M 



where we used (15.11) and (15. 9p . 

From relation (I5.19P we have that J D A s x ip dq = 0. Relation (12. 3 p together withe the above 
estimates imply (I5.20p . The corotational case (I5.2ip also follows since in this case J3 = 0. □ 

H s uniform bounds in the general case. We consider here the general case a(u) = 
Vu. Let us first state the following remark. 

Remark 9. One can easily check that, given four strictly positive constants Ai, A 2 , A 3 , v4 4 we have 
the following property: there exists some u > such that 

A X X 2 + uA 2 Y 2 > A 3 XY + uoA A XY V X, Y 

if and only if A\A 2 > A3A4. Moreover, if the later is true then one can choose lo = 2A ^ A2 ~ A 'i A ± 
which is of the same order as A\A 2 jA\. 

We impose now that the condition above holds true with constants 

OL\ Ot 2 C 2 

Ai = —, A 2 = a 3 , A 3 = Ci— , A 4 = — , 

that is we impose that 

(5.22) cua^R 2 > 2C x C 2 a 2 . 

Let u be as in the previous remark, of the same order as ct\a 3 R A . Assume moreover that 
(5-23) l|A^||| 2( n)+^lKll|^ M <min( I ^,^). 

We multiply (I5.20p by u and add the result to (I5.10p . After using Remark [9] and recalling that 
ij) \ nxD = 4> n we obtain that 

^(||Ai^|| 2 i2(n) + u\\T\\Tr mVL ) + yllA^l^n) + ua 3 \[f 



Ifi 



Let T be the first time such that 



(5.24) ||Ai^(T )||| 2( ^+a;||^(T )||^ =^(^^,^3. 

Then, for t G [0, To], we have that 

l|A^ n (t)|| < " J 



4d 
so that 

d|| A* «»(*) 11^(0) || A^«"(*) ||i, (n) < ^\\A^u n (t)\\l Hn) . 

We also have that 



x m - 2C 2 V u 

so 



1 OtlCts 



wC 2 ||A 2 u ||x3(n)||^ < ||A a u WhIh^ 

< ai \\^v n \\ 2 + ^ 11^11 2 

< -g-||A u Ili»(n) + ""2" I IV Hs^AJ 



We deduce from the above relations that, for £ G [0,T ( 



oj- 



(5.25) IdlAiuii^, + ^m\\- A ) + f|A*«T w + £ o 



which implies that 

||Ai^(T )||i 2(n) +a;||^(T )|||^ < ||A* <||| 2(n) + wKfr^ < min^, ^). 
This contradicts (I5.24p . Therefore the time T cannot exist, so 



l|A^)||| 2(n) +o;||rWII^ M <min( T ^,^) V* > 

and relation (15.251) must hold true for all t > 0. 

We state the result proved in this paragraph in the following proposition. 

Proposition 10. Suppose that cr(u) = V« and that the material coefficients verify relation (I5.22p . 
Moreover assume that 

2 

( 5 -26) IMIir.(n) + C^llwlj^ < min(^, -^). 

Then the sequence u n is uniformly bounded in the space L°° {R + ; H s (Q)) n L 2 (IR+; H S+1 (Q)) and 
the sequence ip n is uniformly bounded in L°°(R+; H^L 2 M ) D L 2 (M+; H^H^). 

Proof. It suffices to show that (15.261) implies (I5.23P for n sufficiently large. This follows at once 
from O, d33D and fl5TT5l) . □ 
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H s uniform bounds in the corotational case. We consider now the corotational 
case a{u) = Vw — (Vu)*. Denoting 



/i = \\A'2 U n \\ L 2 {n) , f 2 = ||A 2 u n \\ L 2 (n) , gi = \\ip Wjjs^ and g 2 = \\ip 
and recalling that ^ \ n D — ' l P n , we have from relations (I5.10P and f)5.2ip that 

and that 

g[ + 2a 3 g 2 < C 2 y / gtf 2 g 2 . 

Using the following two bounds 

C 2 \/g\f 2 g 2 < a 3 g 2 + —^g\f 2 and \fhg~2 < ^-f 2 + ^%02 

4a 3 o z 2 laid* 

we infer that 

(5-27) f[ + ^f 2 < C x hh + gf<7 2 

and that 

(5.28) g' +0l3 g 2 <^L gi f 2 . 

Let £ be a small enough constant to be chosen later but such that 

(5.29) ^(0) < e. 
Let T be the first time such that 

(5.30) 9i(T ) = e. 

We have that giif) < e for all t e [0, T ). Using this in (I5.28|) and integrating in time implies that 
for all t E [0, T ] 

(5.31) gi(t) + a 3 f g 2 <gi(0) + ^ [ f 2 . 

Jo 4a 3 Jo 

Multiplying ( I5.27P by e~ Cl fo h and integrating in time results in 

fi{t) + ^f Q h{s)e C ^>ds < me ^^+^j\ 2 {s)e c ^>ds 
We use now the estimate f)5.3ip above. We infer 



fh < [/.(o) + f^f e^~» fh 

Jo 2a 1 a 3 d 4: 8«ia|d) 4 J 



We now add the following assumption on e: 

(5.32) C!Cl^e eCifrh<ai 

IS 



Assuming that this is true, we further obtain that for all t G [0, T ] 

Going back to (I5.28p . ignoring the second term on the left-hand side and using the Gronwall 
lemma implies now that 

(5.33) 9l (t) <^ 1 (0)exp{^[/ 1 (0) + -^4 I ^ 1 (0)]e^/o^} 

for all t G [0,T ]. If we further assume that 



(5.34) 9l( o)e X p{^[ /l( 0) + ^ 9l (0)]e^} 



< e 



then we observe that setting t = T in ( I5.33|) contradicts (15.301) . We conclude that under the 
hypothesis (I5.29p . f !5.32|) and (I5.34p the time T cannot exist, so all the previous relations hold 
true for all times t > 0. Clearly (I5.29P is implied by (I5.34p . Recalling (15.1 ip we therefore observe 
that there exists some e verifying (15.291) . (I5.32p and (15.341) if we have that 

(5.35) 9l (0)exp{-^L[/ l( 0) + J^m]*^*™*"^*^} 



« 

-da x 7 (KHl^n) + «4||^o|li iL 2 ) 



We state the result proved in this paragraph in the following proposition. 

Proposition 11. Suppose that a{u) = (Vu — Vm)*. There exists a constant C = C(Q,s) such 
that if 

(5.36) I*I^«P{^[MW, + ^WW^^^ 



< n \ ex P 



_ 4 

-Ca x • (||«o||i2 (n) + a 4 \\ip \\ 2 L 2 L2 ) 



then the sequence u n is uniformly bounded in the space L°°(]R + ; H S (Q)) fl L 2 (R + ; H S+1 (Q)) and 
the sequence ip n is uniformly bounded in L°°(R + ; H s x L? hl ) fl L 2 (IR + ; H x H\j). 

Proof. It suffices to show that (15.361) implies (I5.35P for n sufficiently large. This follows at once 
from ([ESD, (J33D and fl5TT5|) . □ 

Remark 12. It is not difficult to see that there exists some constant K = K(Q, s, a x , CK2, CK3, a^, 5) 
such that condition (I5.36P is implied by the following condition: 



^oWmLh < exp 



4 

-K(l + \\u \\ H s {Q) )e Kllu ° r ^ 



End of the proof. Clearly the hypothesis of Proposition [TU] is implied by that of Theorem 
HI and the hypothesis of Proposition [TT] is implied by that of Theorem [2] (see also Remark [TiZl) . 
Therefore, under the hypothesis of Theorem [1] in the general case and under the hypothesis of 
Theorem[2]in the corotational case, we have that the sequence u n is uniformly bounded in the space 
L°°(R + ; H s (tt))nL 2 (R+; H s+1 (ty) and the sequence f n is uniformly bounded in L°°(R + ; H S X L 2 M )C\ 
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L 2 (R + ; H^Hm). Using the equations of u n and f n (13. 4p and (13. 5p this immediately implies some 
time-derivative estimates for u n and f n . A standard compactness argument allows to pass to the 
limit and find a solution (u, f) of such that u G L°°(R + ; H s {Vt)) n L 2 (R + ; H S+1 ({1)) 

and / G L°°(M + ; H^.L 2 M ) fl L 2 (R + ; H^H^). We skip the details since this is very classical and 
straightforward. Next, we clearly have that u is divergence free and tangent to the boundary. 
Passing to the limit in (14. 3[) shows that that relation holds true with /" replaced by /. Finally, 
the equation (13. 5p preserves the sign if the initial data is single-signed so f n > which implies in 
turn that / > 0. This completes the proof of the existence of the solution. 

The uniqueness of solutions is obvious and follows by making energy estimates on the difference 
between two solutions. If (tti,/i) and (w 2 , $2) are two solutions with the same data, then we 
multiply the difference of the equations of U\ and u 2 by Ui —u 2 and the difference of the equations 
of fi and f'2 by 7/ 2 anc ^ ac ^ the two resulting relations. Uniqueness follows easily from the 
Gronwall inequality using that u G L 2 (M + ; Lip) and the linearity in the q variable of the equation 
of /. This is very standard so we skip the details. Remark that even though we show that (12. 2\\ 
holds true in the sense of distributions, i.e. (I2.2p can be multiplied by test functions which are 
compactly supported in the q variable, it can in fact be multiplied by functions which are H]^ in 
the q variable. This follows from a density argument using that C^°(D) is dense in H\ { (D) as was 
proved in [Mas08] . This completes the proofs of Theorems [T] and EJ 

6 Final remarks 

First, we would like to explain here why the condition on the coefficients (I1.5P is necessary in 
the general case on bounded domains. We will observe that a certain cancellation that occurs in 
the case without boundary does not work anymore in the presence of boundaries. When making 
H m estimates on u and H™L 2 M estimates on ip = f — aM we apply d a to the equation of u given in 
(12.11) and multiply by d a u, we apply d a to the equation of ip given in (15.31) (where we dropped the 
superscript n) and multiply by and we add the two resulting relations. We get the following 
right-hand side: 

L l v - ■ L H ■ 9> - h IL v - ■ |va> qM] it 

In the case of a domain without boundary, making an integration by parts implies, after some 
calculations, that the term above vanishes. But in the case of a domain with boundary, the 
boundary terms do not vanish. Moreover, due to the presence of the pressure in the equation of 
u, the first term above should have the Leray projector P in front of V x - making the validity of 
this identity even more unlikely in presence of boundaries. Since the term above does not vanish 
anymore, we need to be able to say that it is small (negligible compared to others) and this in 
turn requires the smallness condition (11. 5ft . 

Second, we would like to explain why the restriction 1 < s < | is necessary. In order to be 
able to control the equation on / we basically need Lipschitz regularity for u. If we assume that 
the initial velocity belongs to H S (Q), then the standard regularity for u obtained through energy 
estimates is L 2 H^ +1 - To get Lipschitz regularity in x we therefore need to assume that s > 1. On 
the other hand, when making the same H s estimates on u n we are led to applying the projector 
P n to the equation of u n (I4.4p and to estimate the right-hand side in H s ~ l . This requires the 
projection P n to be bounded in which implies s — 1 < | so s < |. This explains why the 

condition 1 < s < § is required. We would also like to point out that in dimension three, the first 
requirement that u to be Lipschitz implies s > | while the second requirement does not change 
leading to contradictory assumptions. This means that our approach does not work in dimension 
three. 
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